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Principal Investigator: G. I. Stegeman 
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The enclosed reprints and preprints represent work carried out by personnel sup¬ 
ported under this contract. Summaries of this work along with other work performed 
but not as yet, published are contained within this report. 

We (Stegeman and Seaton) were invited to produce a review paper on Degenerate 
Four-Wave Mixing 1 with guided waves which was published in a special issue of IEEE 
Journal of Quantum Electronics. One of the conclusions of that review was the need 
development and understanding of new nonlinear optical materials suitable for nonlinear 
guided wave devices. 

In light of this need a major part of our research program in the last year has been 
the continuation of research into the semiconductor-doped glass material system. Bulk 
DFWM and photoluminescence 2 experiments on these materials determined that they 
have a large nonlinear index of refraction (n 2 ~ 10“ u m 2 /w) with fast response times 
(r < 16 psec) required for high speed optical signal processing. Along with collaborators 
in the University of Glasgow we successfully developed a successful ion-exchange 
(Na + :K + ) technique for producing low-loss optical waveguides in these materials, the 
details of which are in the enclosed reprint. 2 

We have been studying the effects of incorporating "real" material parameters into 
the operation of nonlinear optical waveguide devices 4 and have determined that the 
most important parameter for guided wave devices is An Mt ie. the maximum change in 
refractive index which can be induced in a material before saturation effects set in. To 
determine this parameter (which has not previously been studied) we have been per¬ 
forming bulk DFWM experiments in a number of semiconductor-doped glass samples. 
By plotting DFWM reflectivity against pump intensity we observe signal saturation over 
a range of samples at ~ 1 MW/cm 2 intensity levels. Using both a two-level atom satu¬ 
ration model and on exponential decay model we determined An sat of 1.7 x 10' 4 and 
1.13 x 10 -4 respectively. Experiments are underway to determine the origin of "per¬ 
manent" gratings which appear to be created via a photo darkening effect in the semi¬ 
conductor glasses. This effect appears to be particularly prevalent in the Schott soda- 
lime glasses which we specially obtained to facilitate the ion-exchange process. The 
problem of "solarization" ie. photo-darkening is known to be particularly prevalent in 
soda-lime based glasses so color filter samples with borosilicate host glasses (ie. BK.7 or 
standard Corning filters) will be further investigated for ion-exchangability. 








. Other work which has been performed with AFOSR support from this contract has 
been of a theoretical nature regarding nonlinear guided wave (NLGW) structures where 
the nonlinearity is sufficiently large to effect the guided wave field distribution. We 
have demonstrated numerically the efficient excitation of stable TEq NLGW's with a 
suitably chosen gaussian input beam. Selective excitation of any of three possible field 
distributions corresponding to the same flux level in a system comprising a linear thin 
film waveguide bounded by two nonlinear, self-focussing media was demonstrated 
numerically. 5 

We also observed that external excitation of a NLGW can produce sequential thres¬ 
hold behaviour through multi-soliton emission from the nonlinear wveguide. 6 This 
effect has intriguing consequences for original optical switching devices. 
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Degenerate Four-Wave Mixing with Guided Waves 
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GEORGE I. STEGEMAN. COLIN T. SEATON, and C. KARAGULEFF 


(Invited Paper) 


Abstract —The theory, experimental implementation, and applica¬ 
tions of degenerate four-wave mixing with guided waves are reviewed. 

I. Introduction 

D EGENERATE four-wave mixing (DFWM) has been 
studied extensively in bulk media for a number of 
years, and many applications have been predicted and 
demonstrated [1]. In the most general interaction geom¬ 
etry. there are three input beams, two of which are coun- 
terpropagating and a third signal beam incident at an ar¬ 
bitrary angle, all of which interact to generate a fourth 
wave via the third-order susceptibility of the medium. The 
process is automatically phase matched and the conjugate 
beam is emitted back along the direction of the incident 
signal beam. The efficiency of the process is proportional 
to the product of the two pump beams, the length of the 
interaction region along the direction of the conjugate 
beam, and the magnitude of the nonlinear susceptibility. 

The electromagnetic waves guided by the interface be¬ 
tween two semi-infinite media, by single or multiple films 
bounded by two semi-infinite media, or by fibers can be 
used to implement DFWM. The key feature introduced 
by guided waves is that the fields decay exponentially 
away from the boundaries into the semi-infinite bounding 
media with Me distances of typically a fraction of the 
wavelength of the radiation being guided. Therefore, the 
effective beam dimension along the direction normal to 
the surface can be on the order of the wavelength of light, 
which corresponds to the minimum beam cross-sectional 
area and hence maximum power density for a given input 
power level. Furthermore, these high fields can be main¬ 
tained for long propagation distances, centimeters in the 
case of thin-film or channel waveguides, and meters or 
kilometers for optical fibers. Thus, the high fields and long 
phase-match distances required for efficient interactions 
can be achieved with modest or low total powers. These 
attractive features were recognized at an early stage and 
have, by now. been applied to an impressive number of 
nonlinear guided wave phenomena [2]. 
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Guided waves, in general, result from a coupling be¬ 
tween an electromagnetic wave and some resonance. The 
resonances can be geometric, for example, as occurs in a 
fiber (or thin-film) waveguide [3], [4] where constructive 
interferences across the core or film result in waveguide 
modes. The resonances can also be material properties re¬ 
lated. For example, the plasma resonance associated with 
the electron gas in a metal is coupled to an electromag¬ 
netic field via the interaction between the field and the 
charges and leads to surface plasmons [5]. Such material 
resonances are usually accompanied by losses and there 
is always a tradeoff between the propagation distance and 
the field confinement [6). For this reason, most work in 
guided wave nonlinear optics for device applications has 
been limited to dielectric waveguides with freely propa¬ 
gating guided waves. 

The first experiments reported [7] on waveguide 
DFWM were performed in CS 2 filled fibers in which the 
long interaction distance was more than enough to com¬ 
pensate for the small nonlinearity associated with CS 2 . 
The phase conjugate properties of the interaction were 
demonstrated. Despite promising applications [8] such as 
wide-angle narrow-band filtering, image-frequency con¬ 
version, and small-sample spectroscopy, there has been 
little progress in phase conjugation in fibers since about 
1980. Recent interest has centered on experiments in 
planar geometries; specifically, experiments using thin- 
film guided waves [9] and surface plasmons [10] were re¬ 
ported. Because of the potential applications to "real- 
time” signal processing [2], [10]. these initial experi¬ 
ments represent the first steps in a field which should de¬ 
velop rapidly over the coming years. 

In this paper, we review progress in degenerate four- 
wave mixing with guided waves. Although the emphasis 
is on experimental results, the basic formulation of the 
problem and its applications will also be discussed. 

II. Theoretical Formulation 

Degenerate four-wave mixing can be carried out with 
either two or four beams guided. Certainly, the most ef¬ 
ficient interactions occur with all four waves guided since 
this configuration optimizes the power density in each 
wave for a given power. For the two-wave guided case, 
two externally incident waves can be used to phase con¬ 
jugate a guided wave [12]. or two guided waves can be 
used to phase conjugate an externally incident plane wave I 
field [13]. Here we consider only the four guided wave 
case. 
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A. Guided Waves 

Guided waves are electromagnetic fields which satisfy 
both Maxwell’s equations in every medium into which the 
fields penetrate and continuity of the tangential boundary 
conditions at every interface. For the simplest case of iso¬ 
tropic waveguide media with propagation along the x axis, 
the sth input wave (s = 1 , 2, 3) in the mth normal mode 
an be written as 


£(r, t) = 4 E ims \r) a im s '(x) 


■ exp [t(wt - 0 {ms % ■ r)] + c.c. ( 1 ) 

where & m s) is the effective index obtained from the so¬ 
lution of the appropriate dispersion relation and E m , \r) 
is the corresponding normal mode field. Since these so¬ 
lutions are normal modes, they satisfy the orthogonality 
relations 



dr, E {m \r) • £«-(r) = 5 m .„ 


( 2 ) 


where the vector r f is orthogonal to the propagation direc¬ 
tion defined by k and A m n is the normalization constant 
chosen so that |a ( " I) (x )| 2 is the guided wave power for two- 
dimensional waveguides or the guided wave power per 
meter of wavefront for planar guided waves. These 
expressions are valid for all guided wave fields, fibers, 
channel waveguides, surface plasmon waves, etc. 

The mixing of three input guided waves leads to non- 
near polarization fields which are usually expanded in 
products of the mixing fields. A nonlinear polarization 
field of frequency oi is generated in each of the media in 
which the guided wave fields exist. It is usually written 
as 


P( N r L „ = I exp [i(wr - 0 p k p x)\ + c.c. (3) 


where 0 p k p is the polarization wavevector. The generation 
of the nonlinearly produced conjugate wave can be treated 
' • the slowly varying amplitude approximation using well- 
lown guided wave coupled mode theory [3], [4] with the 
result that, for propagation along the x axis. 



(x) = i 


W|XpC 

20 in!, A z „ n 



dr, P% ■ E* {n s \r) 


• exp [-i(0 p - 0 ( "' J, )x). (4) 

For a fiber or a channel waveguide, the integral is taken 
over the cross-sectional area, and in the case of planar 
waveguides, the integral is performed ove. the depth co- 
"dinate. 



Fig. 1. The planar four-wave mixing geometry discussed in this paper. 
The pump beams are I and 2. the signal wave is 3. and the conjugate 
wave is 4 [14). 


with wavevector [0 < " ,4, *4] as shown in Fig. 1 for a planar 
guided wave geometry. Wavevector conservation is given 
by 

IS*"’ 0 *, + & m '- 2) k 2 + & m ' h k } + = 0 (5) 

with *| = = * 4 . If the two input beams 1 and 2 

are exactly contradirectional and have the same mode 
number m = m\ 0 (ml> *, = -& m - 2, k 2 . Therefore, it is 
necessary that m" = n so that ^ nA) k t = -0' n - 3) k i , and 
signal beam 4 is generated backwards along the direction 
of incidence of beam 3. independent of the incidence di¬ 
rection of beam 3. (For a fiber, all four beams propagate 
along the fiber axis.) Note that this process does require 
that beams 1 and 2 have the same mode order m, that 
beam 4 is radiated into the same mode index n as beam 
3, and m and n are not necessarily equal [14], Because 
there are multiple guided wave solutions associated with 
the same guided wave polarization characteristics, the 
spectrum of possible interactions is larger than that for 
plane waves. (The other consequence is that the possible 
directions of the incident waves are restricted to one (fiber 
or channel waveguide) or two (planar waveguide) dimen¬ 
sions for an all-guided wave interaction.) 

For an isotropic material, the conjugate polarization 
term [14] is given by 

P ^*(r) = [4<o\ °22( — w, u>. ~w. u>) 

• £*'"V> Ej m :, (r) E) mA \r ) 

+ 2toXi::i( —w. «. -w. u) E\ m "(r) 

■ Ej m 2> (r) £*'"• V)) a,(x) a : (x) af(x ) (6) 

where the notation for the field amplitudes has been sim¬ 
plified so that a^x) ■ a'" "(x). This expression simplifies 
far from any resonant behavior in the x l3) terms and 
Xna = X i 22 i = Xeff- For guided wave fields. 

P NL ”(r) = 2«5n 2 n :> [$ £?,'"'V) ££ V) ££- J '(r) 

+ \ E?; n 'V) £;7 : ’(r) ££ "(r» 


B Degenerate Four-Wave Mixing in Guided Wave 
Geometries 

The same formalism can be used to discuss the four- 
wave mixing in one- and two-dimensional waveguides [ 8 ], 
(14]. There are three input beams with wavevector 
[ 0 °""*,], [ 0 *" 2) * 2 ], and ( 0 *'" ^itj] and one output beam 


• a,(x) a : (x) at(x) ( 7 ) 

in each medium (labeled 7 ) where n :> * 0. Here we have 
defined 3xefr = where n is the intensity-de¬ 

pendent refractive index for the yth medium. Clearly, 
similar equations can be derived for the nonlinear polar¬ 
izations driving the pump and signal beams. 


i 

1 
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Here we consider only the two cases which have been 
studied experimentally. For fibers or channel waveguides, 
the coupled mode equations are 

~~ a 4 (x) = f( a )1 22 + * 12 : 1 ) a*(x) G'S*) d|U> (8a) 
ax 

d 

— a } (x) * —((a 1122 + * 1221 ) a*(jO a\(x) a 2 (*) (8b) 

dx 


with 


*1122 — 


2u>to 

3*[jS 


J dr, n;n 2y [E*l" 4 \r) £*'"- 3, (r> 


r (w -2)/ 


*1221 


V) £;?-»! 

= J ^ dr, n;n 2y [£*'"V) ^‘"V) 


(9a) 


£y‘ 


(m. 1) 


(r) £yj ,m - :, (r)J. 


(9b) 


In the weakly guiding case, that is, one m which the index 
difference across the core-cladding interface for fibers or 
across the channel-bounding media interfaces is small. 
a 1221 - A) ,22 and (8) simplify to the well-known equations 
for four-wave mixing with plane waves with, of course, 
a 1122 given by the •‘overlap integrals.” In the absence of 
pump depletion with the boundary conditions u 3 (0) = a, 
and a 4 (L) = 0 where L is the length of the interaction 
region, the results are 


a 3 (x) = a, , tan — sin — + cos —; 

L *cr cr *crJ 


(10a» 


a 4 (x) = ia* 


L x 
tan — cos — 



— = (a,, 22 + *1221) a i( 0 ) a 2 ( 0 ). 

*rr 


(10b) 

( 11 ) 


Consider also the case of waves guided by planar sur¬ 
faces. Assuming that TE beams 1 and 2 propagate along 
x ’ which is oriented at an angle 0 to the x axis along which 
TE beams 3 and 4 are traveling, 

£■<«. ior 2 )( r ,> m E' n n (r’) [sin 0, cos 0. OJ; 
E'” im *\r) = £"°Yr) [0. 1, OJ. (12) 


Again, in the limit of no attenuation or depletion for the 
input beams, and assuming that beams 1. 2. and 3 are 
introduced into the interaction region at x ' = 0, x‘ = L ‘. 
and x - L, respectively, substituting into (8), and inte¬ 
grating the signal beam from 0 — l gives the same form 
as (10) and (11), but with 


*1122 + *1221 


3/i;.„r' 


(2 ■+■ cos 2 0) 



dz n y (z) 2 n 2y (z) 


■ [£';•" <:) £',"■ 2) 


(z) Ef in i \z)E*'" *'(z)]. 


(13) 


In the general case, it is necessary to take incident beam 
depletion into account which leads to a series of coupled 
mode equations between the amplitudes of the various 
beams. 


III. Experiments in Fibers 

In 1975, Yariv [15] suggested using nonlinear mixing 
as a means of compensating for phase abberations in¬ 
curred in transmitting an image down an optical fiber. He 
subsequently pointed out the possibility of performing this 
phase conjugation by degenerate four-wave mixing in a 
multimode waveguide [16] and suggested using a hollow 
core fiber filled with CS 3 as the nonlinear mixing medium. 
Experiments based on this idea were reviewed by Hell- 
warth in 1982 [17], and to our knowledge. there have been 
no further results published since then. Thus, the only ex¬ 
periments to have been done are those of Jensen and Hell- 
warth [7] and Au Yueng er al. [18] in CS 2 filled fibers. 
The first authors achieved high-fidelity conjugation with 
multimode pump waves using a Q-switched ruby laser and 
achieved DFWM reflection efficiencies of =50 percent 
for 35 kW peak pump beam powers. They also demon¬ 
strated that the wave being conjugated could serve simul¬ 
taneously as its own pump wave. Au Yueng et al. [18] 
used single-mode pump waves to obtain 0.45 percent 
DFWM conversion efficiency for only 6 mW of CW pump 
power from an argon ion laser. However, the phase con¬ 
jugate nature of the DFWM signal was not ascertained. 

In related experiments, image transmission through a 
multimode optical fiber has been demonstrated [19] by 
using a photorefractive BaTiO. crystal to correct for modal 
distortions external to the fiber. Stimulated Raman [19] 
and Brillouin [20] scattering in fibers can also be used to 
achieve phase conjugation, but these topics are outside 
the scope of this paper. 

IV. Experiments with Planar Glided Waves 

Two experiments on four-wave mixing in integrated 
optics waveguides have been reported to data [9], [21]. 
The first deals with optically induced damage in lithium 
niobate waveguides [21]. A TE guided wave was created 
due to stray scattering via the photorefractive effect, and 
the mixing of the incident TM wave with this TE light 
creates a grating which further enhances the TM — TE 
conversion. Since the gratings formed are written more or 
less permanently into the waveguide, this does not cor¬ 
respond to the usual four-wave mixing phenomenon. The 
characteristic time scale was typically seconds. 

The second experiment was performed using freely 
propagating guided waves [9]. The nonlinear medium 
CS 3 . which was used as a cladding on top of a thin-film 
waveguide, was located inside a cell optically contacted 
to the surface of a thin-film glass waveguide, as shown in 
Fig. 2. Three coupling pnsms were required for the four 
beams needed in the interaction, leading to very difficult 
alignment problems. Since the nonlinear mixing occurred 
via the evanescent tails of the three incident guided waves 
and because the nonlinearity of CS : is small (n 2l . = 3 x 
10''* m 2 /W). the reflectivity (fractional conversion of 
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Fig. 2. Sample and coupling geometry for degenerate four-wave mixing 
in a thin-film waveguide showing the (hin-film waveguide (Coming 7059 
sputtered onto a glass microscope slide), the arrangement of high-index 
coupling prisms, and the glass cell for holding the CS,. The probe beam 
intersects the pump beams at 90° under the liquid cell (9). 


Fig. 3. The application of four-wave mixing to (a) convolution (of wave¬ 
forms 1 and 2) and (b) time inversion lof pulse 1) [II]. 


input pulu 
* at k >0 532 



beam 3 into beam 4) was only = 10 -9 . Nevertheless, the 
cubic dependence of the four-wave mixing signal was ver¬ 
ified experimentally for 15 ns long pulses from a Q- 
s witched, frequency-doubled Nd: YAG laser. The ob¬ 
served signal was consistent with calculations based on 
the reorientational nonlinear mechanism. 

V. Planar Waveguide Applications 

This interaction has some possible applications to real¬ 
time all-optical signal processing [11], two of which, con¬ 
volution and time inversion, are illustrated in Fig. 3. The 
two input signals to be convolved are beams 1 and 2 which 
have pulse envelopes C/,(f - x/v) and (A(f + x/v). 
Whether the convolution process occurs or not can be 
controlled by beam 3 which is assumed to have a constant 
amplitude during the overlap of pulses 1 and 2 inside the 
waveguide. Assuming that 6 - t/2 for this case (geom¬ 
etry in Fig. 1) and that L is small relative to the charac- 
eristic pulse size (to preserve detail in the convolved sig¬ 
nal), the fourth beam signal radiated is directly pro¬ 
portional to the instantaneous overlap of beams 1 and 2. 
Mathematically, the time evolution of the signal beam can 
be written as U t (t) «{£/,([- 2r) tA(7) U^dr where it 
has been explicitly assumed that the control beam is 
"turned on.” For this application, the input signal length 
cannot be larger than the waveguide dimension which 
limits the process to pulses of 100 ps or less. Note also 
the time compression by a factor of two. 

Another potential application of four-wave mixing to 
real-time processing is time reversal of an optical wave¬ 
form, as illustrated in Fig. 3(b). Here 9 = 0 and beams 
1 and 2 are the waveform to be reversed and a control 
beam of constant amplitude, respectively. Beam 3 is a 5 
function in time (very short pulse) and 4 is the time re¬ 
versed signal. A crucial factor is that beams 1 and 3 can 
be propagated in different modes. If the mode number of 
beam 3 is larger than that of beam 1, beam 3 can overtake 
beam 1 from "behind.” Thus, beam 4 is "read out” from 



Fig. 4. Experimental setup used for the observation of optical phase con¬ 
jugation with surface plasmons. The inset shows (he waveveciors of the 
three incoming waves in the plane of the interface (BS beam splitters. 
VDL: variable delay lines. PA: pnsm-silver film assembly mounted on 
a rotary table. L: 1 m focal length lens. A: 0 5 mm diameter aperture. 
PMT: photomultiplier tube. PD photodiode (101 l 

trailing edge to leading edge, which corresponds to time 
reversal of the pulse envelope. 

VI. Surface Plasmon Experiments 
A variety of experiments have been earned out by Nunzi 
and Ricard [10] on phase conjugation with surface plas¬ 
mons in attenuated total reflection (ATR) geometries. In 
the visible and near infrared, the propagation distance for 
surface plasmons is typically tens of micrometers and the 
nonlinear mixing takes place under a single coupling 
prism. 

In the first set of experiments, three incident beams were 
coupled into surface plasmon modes as shown in Fig. 4 
and the conjugate beam was sampled via a beam splitter. 
For a silver film on a glass prism probed in the Krctsch- 
mann geometry, the four-wave mixing signal observed 
was due to heating of the metal film. This was verified by 
pump-probe experiments on the large spacing nonlinearly 
produced grating. With pulses 28 ps long, 0.6 ns decay 
times were measured, in good agreement with calcula¬ 
tions. When liquid ethanol was placed between the cou¬ 
pling prism and the silver film, thermal gratings were 
formed in both the liquid and metal, and interference ef¬ 
fects were observed. The conjugate reflectivity in these 
experiments was on the order of 10 -4 . This case has been 
analyzed not only by the experimental group [10], but also 
by Ujihara for nonlinearity in the dielectric medium [22] 
and in metal [23] (modeled as a free electron gas). 
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• In the second set of experiments, two counterpropagat- 
ing surface ptasmons were excited under a coupling prism 
to produce a refractive index grating. This grating was 
then used to couple in a third wave, incident from the air 
side of the metal film, into a surface plasmon mode. The 
nonlinearly in-coupled surface plasmon was observed 
when it was coupled out by the prism. The mixing of the 
two counterpropagating surface plasmon modes was also 
found to generate high-frequency acoustic waves. 

■ There has been some speculation that the efficiency of 
DWFM by surface plasmon could be increased by using 
the long-range mode guided by metal films of thickness 
on the order of the optical field skin depth (6], (14). Al¬ 
though it appears theoretically that the efficiencies should 
be larger (6), (14], it has been argued by Nunzi and 
Ricard (10) that using such thin films would not eliminate 
the heating effects which dominated their experiment. 

VII. Summary 

The single factor which has held ba^ k the development 
of degenerate four-wave mixing in guided wave structures 
has been the lack of nonlinear materials suitable for fab¬ 
rication in waveguide form. In the last few years, a num¬ 
ber of new material systems have been identified which 
can be fabricated into thin film, and probably fiber form. 
They are essentially semiconductor-(24]-[27] and metal- 
doped (28) glasses and thin organic films (29). All of these 
materials have very fast relaxation times, with the semi¬ 
conductor glasses being the slowest with response times 
(electronic) of tens of picoseconds. At the present, the 
semiconductor-doped glasses are the only ones with suf¬ 
ficient nonlinearity to result in waveguide four-wave mix¬ 
ing with subwatt peak power levels. However, improve¬ 
ments are expected in the field of nonlinear organic 
materials so that they become competitive with the doped 
glass materials. 

With such highly nonlinear and ultrafast materials, 
some of the promising applications to all-optical signal 
processing should become feasible, and we look forward 
to these possibilities with great anticipation. 
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T. J. Cullen and C. N. Ironside 

Department of Electronics and Electrical Engineering, University of Glasgow, Glasgow, G12 8QQ, 

United Kingdom 

C. T. Seaton and G. I. Stegeman 

Optical Sciences Center. University of Arizona, Tucson, Arizona 85721 

(Received 13 August 1986; accepted for publication 29 September 1986) 

The properties of planar waveguides and directional couplers fabricated by potassium/sodium 
ion exchange in semiconductor-doped glasses are described. 


A new spectrum of ultrafast (picosecond) all-optical 
devices becomes possible when an optical waveguide con¬ 
tains a medium characterized by an intensity-dependent re¬ 
fractive index. 1 Although many interesting device concepts 
have been proposed, 1 progress has been hampered by the 
lack of nonlinear materials with both large nonlinearities 
and picosecond response times, and from which optical 
waveguides can easily be fabricated. To date, only a nonlin¬ 
ear coherent coupler, first proposed by Jensen, 2 has been 
demonstrated with cw excitation in strain-induced wave¬ 
guides in GaAs/AlGaAs multiple quantum well material. 3 ' 4 
This device consisted of two parallel, channel waveguides 
and, when light is initially inserted into one channel, the 
transfer efficiency into the second channel depends upon the 
initial power in the incidence channel. A further refinement 
is a nonlinear Mach-Zehnder interferometer as described by 
Lattes et alS Recently it has been shown that semiconduc¬ 
tor-doped glasses (color filters), which are glasses contain¬ 
ing small (10-100 nm) crystallites of CdS^Se, _ x .have use¬ 
fully large nonlinearities 6,7 (n 2 = 10“ 14 m 2 /W) for guided 
wave devices, coupled with picosecond relaxation times 7 
(r<10 _ " s). For example, the critical power P e foranon- 
linear coherent coupler 2 acting as an all-optical switch 
should be approximately 1 W for a channel cross-sectional 
area of 10 /xm 2 and a linear transfer length of 3 mm. Since 
such glasses are relatively inexpensive and can be made with 
10% concentrations of sodium, they are promising candi¬ 
dates for waveguide fabrication by ion exchange, one of the 
more simple waveguide fabrication technologies.*' 9 In this 
letter we describe the fabrication and the waveguiding char¬ 
acteristics of both planar waveguides 10 and directional cou¬ 
plers in semiconductor-doped glasses. 

The color filter glasses used in this study consist of a 
bulk soda-potasb-lime glass matrix doped with small crys¬ 
tallites of the semiconductor compound CdS z Se, _ z with x 
chosen to give a red colored glass with a band edge around 
590 nm. The nonlinearity is resonant with the band gap and 
this wavelength is convenient for subsequent characteriza¬ 
tion of the nonlinear properties of the devices with a mode- 
locked rhodamine 6G dye laser. The glass content is 9.5% by 
weight sodium and 3.5% by weight potassium. The sub¬ 
strate index was measured using an Abbe refractometer to be 
1.5185 ± 0.0001 at 632.8 nm. Planar waveguides were fabri¬ 
cated by immersing the filter into a bath of molten potassium 
nitrate at 340 *C (this temperature just above the melting 
point of 334 ‘C was found to give the most reproducible re¬ 


sults ) in which sodium ions leave the glass to be replaced by 
potassium ions from the melt. This exchange results in a 
region of higher refractive index An of around 0.009 at the 
surface which decreases monotonically with depth. To fabri¬ 
cate stripe waveguides, an aluminium diffusion mask is de¬ 
posited on the glass: the mask has apertures through which 
selective ion exchange can take place thus defining the wave¬ 
guide dimensions. 

The fabrication of single-mode planar and channel 
waveguides requires a detailed investigation of the ion-ex¬ 
change process. A characterization study of K*"/Na + ion 
exchange in color filter glass was undertaken similar to those 
reported in Ref. 8 (Ag + /Na + ion exchange) and Ref. 11 
(K*'/Na + ion exchange) for soda-lime microscope slide 
glass and Ref. 12 (K^/Na + ion exchange) for BK-7 glass. 
A number of planar waveguide samples were prepared in 
which the diffusion times t were varied between the range 15 
min to 24 h. For every sample the effective index of each 
mode was determined by the prism coupling method to an 
accuracy of ± 0.0001 with a He-Ne laser operating at wave¬ 
length of 632.8 nm. The results for both TE and TM polari¬ 
zations are shown in Fig. 1. 

The theoretical dispersion curves which produced an 
excellent fit to the data were derived as follows. First, the 
inverse Wentzei-Kxamers-Brillouin (WKB) method 13 was 
used to obtain the mode turning points and also the surface 
refractive index n,, consistently found to be 1.5272 for TE 
modes and 1.5305 for TM modes. In order to find depth 
index profiles which accurately fit these points, the finite 
difference method was used to calculate solutions to the non¬ 
linear diffusion equation in one dimension given by 


dc _ d f 

D , dc\ 

(la) 

dt dx \ 

1 —ac dx)' 

C = C./Co, 

a - 1 - (Z>,/2> 0 ) , 

(lb) 


where c is the normalized potassium ion concentration (c, 
being the concentration of the potassium ions and c 0 the 
total number of diffusing potassium and sodium ions in the 
glass) and D , and D 0 are the self-diffusion constants of the 
potassium and sodium ions, respectively. The surface con¬ 
centration ratio was determined by laser induced ion mass 
analysis (LIMA) of the exchanged glass surface layer to be 
9:1 K + /Na*, giving a boundary condition of c = 0.9. The 
refractive index is assumed to vary linearly with the potas¬ 
sium ion concentration. The shape and depth of the index 
profiles are adjusted by varying the parameters a and D ,: 
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FIG. I. Theoretical dispersion curves calculated using the Wentzel- 
Kramers-Bnllouin (WKB) method for the optimum index profiles com¬ 
pared to measured mode effective index values for samples fabricated at 
T =* 340 ‘C. (a) TE modes and (b) TM modes. 


values of a = 0.3 and = 4.95x10*“ m 2 s~‘ were 
found to give the best least-squares fit for the TE modes, 
while a = 0.3 and D , = 4.50X 10 ~ 14 m 2 s ~ 1 give the best 
fit for the TM modes. The apparent anisotropy of ion-ex¬ 
changed waveguides, as suggested by the discrepancy 
between the diffusion constants and the surface index values 
for the TE ana TM modes, is thought to be due to stresses 
introduced into the glass by the different sizes of the two 
ionic species being exchanged. 13 

Using the best fit refractive index profiles and defining 
d , to be the depth where the index change is equal to 1/e of 
the surface index change, a linear relationship between d, 
and t 1/2 was obtained. From the slope, the “effective” diffu¬ 
sion constant [defined by the relation (d, = {D f t) w2 ] is 
given by D, = 10.80X 10" 16 m 2 s~‘ for TE modes and 
D, =9.32x10* 16 m 2 s* 1 for TM modes. Finally, the 
WKB method was employed using the diffusion parameters 
obtained above to calculate the curves fitted to the data in 
Fig. 1 . The excellent fit indicates that the refractive index 
profiles and the D t values determined are reasonably accu¬ 
rate. This agreement is particularly important for predicting 
the short diffusion times required for single-mode wave¬ 
guides. A diffusion time of 1 h was chosen from these results 
as the optimum time for the fabrication of single-mode wave¬ 
guides. 
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FIG. 2. Schematic diagram of directional coupler pattern used to fabricate 
devices. This pattern was repeated three times on the mask for different 
waveguide separations (d) of (a) 4.0, (b) 4.5. and (c) 5.0/am. The interac¬ 
tion length was varied in steps of 0.5 mm over the ranges of (a) 1.0-3.00 
mm, (b) 2.0-4.0 mm, and (c) 3.0-5.0 mm. The waveguide width (w) was 
nominally 4.0 ftm for all the couplers on the mask. 


The theoretical effective index values of stripe wave¬ 
guides and the transfer length for directional couplers were 
calculated using the numerical model derived by Walker 
et al First, the model is used to solve the nonlinear diffu¬ 
sion equation in two dimensions to give the refractive index 
profile n(x,y), then the Helmholtz wave equation is solved 
for such a region using the variational metnod (the model is 
restricted to quasi-TE modes) to give the mode effective in¬ 
dex values. 

For a diffusion time of 1 h at 340 *C the slab effective 
index is 1.5205. The theory predicts that stripe waveguides 
formed with the same diffusion parameters will be single 
mode for aperture widths w of 3-5 /im with effective index 
values from 1.5189 to 1.5196. The mask design used to ex¬ 
perimentally verify directional coupler calculations is shown 
in Fig. 2. 

Several samples were fabricated using this mask and the 
edges polished to facilitate end-fire coupling (this process 
removed 0.5 mm from the interaction length of each device). 
Light from a He-Ne laser was coupled in turn into each de- 
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INTERACTION LENGTH 
Lj (mm) 

FIG. 3. Plot of the measured normalized coupled power < I — P. P ,) vs 
interaction length (£,) for the three sets of similar couplers on a sample. 
The best St cos 2 1 irz/2L c ) curves are also shown. The results are shown in 
Table I. (The lengths are 0.5 mm less than the lengths in Fig. 2 because this 
amount is lost in the polishing process.) 
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TABLE 1. Comparison of theoretical and measured coupler transfer length 
for devices prepared at r = 340 ‘C for 1 h. 


Coupler 

set 

Waveguide 
width u ■: 
(pm) 

Waveguide 
separation d 
(M"l) 

Transfer length L, (mm) 
measured theoretical 

a 

4.2 

4.2 

3.5 ±0.25 

3.5 

b 

4.2 

4.6 

5.2 ± 0.25 

5.5 

c 

4.2 

5.3 

7.7-0,35 

7.9 


All guide widths and separations measured to an accuracy of ±0.1 fim 


vice on a sample at the single-guide end and the power in 
each of the output guides was determined by focusing the 
near-field image onto the camera of a calibrated vidicon sys¬ 
tem. Figure 3 shows the power remaining in the input guide 
versus the device interaction length for the three sets of de¬ 
vices. The best fit cos 2 ( rrz/2L c ) curves are also shown to 
estimate the coupling length. Table I gives the actual device 
parameters and measured transfer lengths compared with 
the predicted transfer lengths calculated from the theoretical 
model. Very good agreement was obtained. 

In summary, good quality single-mode directional cou¬ 
plers can be made by K * /Na * ion exchange in semiconduc¬ 
tor-doped glasses. The waveguiding properties of a direc¬ 
tional coupler can also be predicted with the model 
presented here. Experiments are currently under way to 
measure the nonlinear properties of the directional couplers, 
as well as the Mach-Zehnder devices. 
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Gaussian beam excitation of TEo nonlinear guided waves 
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Selective, efficient excitation of TEo nonlinear guided waves is demonstrated numerically. For 
a thin film bounded by two self-focusing media, three different field distributions 
corresponding to the same flux level can be excited independently by suitable Gaussian input 
beams. 


To date several integrated optics devices have been pro¬ 
posed utilizing nonlinear waveguides (NLWG’s), for exam¬ 
ple. upper and lower threshold devices, nonlinear couplers, 
optical switches, etc. 1 Such proposals have invariably been 
based on the NLGW dispersion curves (guided wave flux 
versus effective index) 2 even in the absence of data on the 
stability of the various branches and of the possibility of ex¬ 
ternally exciting the NLGW’s. Moloney et aI. 3A and Jones 
and Moloney 3 have recently addressed the stability question 
both numerically and analytically, and a comprehensive pic¬ 
ture is now emerging. In contrast, the problem of external 
excitation of NLGW’s (in the waveguide context) has not 
been discussed so far, and is the subject of this letter. Specifi¬ 
cally, we show that, for a fixed input flux, selective excitation 
of the corresponding TE^ NLGW’s by end firing of appro¬ 
priately tailored Gaussian beams is indeed possible. Implica¬ 
tions for device design concepts are also discussed. 

We treat the case of TE waves of frequency to in a planar 
NLGW as shown in Fig. 1. For this preliminary study the 
electric field is assumed homogeneous in the y direction. 
Then writing the electric field as 

E(r.t) = \y[W(xj)e“ *—*’+ c.c.] (1) 

gives in the usual slowly varying envelope approximation, 

HP— + 10 2 -n 2 (x,\fVI 2 )]r = 0, (2) 

02 ax' 

where the refractive index in the various media is described 
by 

/i J (x,|#'| J ) = /t 2 r +a Y \W(xj)\\ r = c,/4, (3) 

and the subscripts c, /, s refer to the cladding (jc< — d), 
film ( — d<,x<,d), and substrate (x>d), respectively. Here 
the coordinates x and z, and d are in units of k 0 ~ 1 = c/co. The 
NLGW’s of the system are solutions of Eq. (2) such that 
W(xs) = fV 0 (x). Equation (2) along with the boundary 
conditions of continuity of W and dW /dx at the interfaces 
(jc = ±d) can then be used to generate the dispersion 
curves for the system. Note that due to the nonlinearity as¬ 
sumed in this problem, several NLGW solutions may exist 
for a given 13 1,2 An example of the dispersion curves for the 
TE,, NLGW’s of a waveguide with nonlinear cladding and 
substrate, and a linear film, is shown in Fig. 2 (we use the 
parameter values given in this figure throughout). We 
choose self-focusing nonlineanties since this case contains 
by far the richest set of new phenomena. 12 The regions 
marked by the dashed lines are unstable under propagation. 


and the remaining regions are stable and are labeled I, II, and 
III, respectively. 5 The inset beside each stable region indi¬ 
cates the field distribution of the NLGW on that branch. For 
example, in region II the NLGW is associated with self-fo¬ 
cusing in the medium with the lower nonlinear coefficient, 
whereas in III it is associated mainly with self-focusing in the 
other bounding medium. At high flux levels these solutions 
degenerate into nonlinear surface polaritons at their respec¬ 
tive interfaces, each independent of the other interface. 

To investigate the excitation of NLGW’s we have solved 
Eq. (2) with initial data corresponding to a collimated 
Gaussian input beam; 

= (4) 

where W N is used to adjust the input beam flux. x Q is the 
displacement of the beam center from the guide center, and 
w 0 is the Gaussian spot size. This boundary condition corre¬ 
sponds to end firing a Gaussian beam onto the NLGW. Note 
that we ignore any nonlinear reflection that may occur at the 
air-waveguide interface. This is valid since the input beam is 
incident normally onto the interface whereas nonlinear ef¬ 
fects occur only for angles close to the critical angle for total 
internal reflection. 5 7 

For our numerical experiment we have fixed the input 
flux (S, n = 0.165). and as indicated by the horizontal line in 
Fig 2, there are three distinct, stable NLGW solutions. For 
each solution we examined the NLGW profile and estimated 
values for x 0 and w 0 such that the input Gaussian would 
resemble the corresponding NLGW. For example, for the 



FIG I Nonlinear waveguide geometry, a linear thin film of thickness Id is 
bounded by two self-focusing nonlinear media 
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FIG. 2. Guided wave flux ( S) vs effective index ( P) for two seif-focusing 
bounding media. The solid and dashed lines indicate stable and unstable 
regions, respectively, and the insets show the nature of the NLGW’s on the 
various branches. Parameter values are n c — n, =» 1.S5, n f m 1.37, 
a, » 10- 1 . a, - 8x 10-’, a, - 0, andd - g. 

solution on branch I of Fig. 2 we chose iv 0 = 10.0 and 
x„ = 0.0, whereas for branch II we have w 0 = 6.0 and 
x 0 = 10.3. The evolution of the Gaussians for branches I—III 
are shown in Figs. 3(a)-3(c), respectively. In all cases the 
Gaussians evolve toward the NLGW they were tailored to 
excite. This was verified by detailed examination of the final 
beam profile and the corresponding NLGW. We have there¬ 
fore demonstrated that selective excitation of the TEo 
NLGW’s at a fixed flux level using a Gaussian beam is possi¬ 
ble. These results indirectly imply that selective excitation is 
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FIG. 3. Evolution of Gaussian input beams for S m > 0.163 and (a) x„ = 0, 
u>„ m I0.O. (b) x 0 « 10.3, ui„ — 6.0; (c) x„ x —10.6, ui 0 » 4.0; (d) 
x, » - 9.0, u/ 0 a 4.0. The propagation coordinate is in units of free-space 
wavelengths, and the vertical lines are the waveguide boundaries. 


also possible if only the cladding is nonlinear, in that case the 
only difference is that branch II vanishes and there are only 
two NLGW’s. We have furthermore verified numerically 
that Gaussian beam excitation is possible for many values of 
the ratio a c /a,. The results presented here are representa¬ 
tive of our general findings. 

A new question immediately arises from these results, 
namely, how sensitive is the coupling to the initial Gaussian 
beam profile? First, we stress that our estimates for z 0 and u> 0 
came from visual examination of the NLGW's which indi¬ 
cates that there is no exceptional sensitivity to initial condi¬ 
tions. However, if the input beam profile is too different from 
a NLGW it decays mainly via emission of radiation and/or a 
soliton which is ejected into the nonlinear medium with the 
larger nonlinear coefficient. 4 (If both media are identical, 
two solitons are emitted, one into each medium.) An exam¬ 
ple of the soliton decay is shown in Fig. 3(d) which corre¬ 
sponds to the situation in Fig. 3(c) but with the beam center 
misaligned from x 0 = — 10.6 tox 0 = — 9.0. More details of 
these decay processes will be published elsewhere. 

These results show that the input beam profile is a signif¬ 
icant control parameter in considerations of NLGW excita¬ 
tion. The input flux alone is not sufficient to characterize the 
input This has implications for device design concepts. For 
example, previous proposals for optical power limiting using 
NLGW's have relied on the fact that for certain configura¬ 
tions (e.g., a waveguide with a self-defocusing cladding 
which starts above cutoff), there exist no NLGW’s above a 
certain critical flux. 1 For the case shown in Fig. 2, we can 
couple a Gaussian beam (x 0 = 0.0, w 0 = 10.0) into branch I 
very efficiently. (The system is essentially linear on this 
branch.) However, for input fluxes above which this branch 
does not exist, our input beam does not match the solutions 
on the two remaining branches (II and III) and the input 
decays via emission of radiation and/or solitons, leaving a 
NLGW on branch I of reduced flux. This process produces 
threshold behavior and optical limiting. Without proper 
consideration of the input profile, this action cannot be un¬ 
derstood since NLGW’s exist for all possible input fluxes. 

In conclusion, we have demonstrated that selective exci¬ 
tation of stable TEo NLGW’s using a Gaussian input beam is 
possible. Indeed it seems that simply matching the input 
beam to a particular NLGW is sufficient to guarantee selec¬ 
tive excitation. 
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We demonstrate numerically that external excitation of a nonlinear waveguide can produce 
sequential threshold behavior via multisoliton emission from the waveguide. This behavior is simi¬ 
lar to that predicted to occur at a nonlinear interface. 


Following the pioneering work of Kaplan 1 on the re¬ 
flection of a plane wave from a linear-nonlinear interface, 
several papers appeared on this subject. 2 ~ q This interest 
stemmed from the possibility of bistable reflection from 
such an interface. 1 Although it is now generally accepted 
that bistable reflection is not possible (as opposed to a 
hysteretic response) several of the results from the theoret¬ 
ical studies, such as the nonlinear Goos-Hiinchen 7 q effect 
and the transmission of self-focused channels (or solitons) 
through the interface, 2 7 are of interest in their own right. 
To the best of our knowledge only the former effect has 
received further attention. 

In this Brief Report we consider the external excitation 
of a nonlinear waveguide (NLWG) where the Film and 
substrate are linear but the cladding displays a nonlinear 
refractive index (optical Kerr effect). Intuitively one may 
expect that effects reminiscent of the nonlinear interface 
problem can arise since the NLWG comprises at least one 
such interface. Here we report results showing the 
transmission of solitons through the Film-cladding inter¬ 
face, the number and angle of emission of which depends 
on the input flux. [The input beam proFtle is held Fixed as 
the zeroth-order transverse-electric (TE) mode of the 
linear guide ] The flux trapped in the waveguide then 
shows sequential threshold behavior as a function of the 
input flux as a result of this multisoliton emission. An in¬ 
tuitive explanation for the thresholds utilizing the NLWG 
dispersion curve for the system is given. We remark that 
the emission of a soliton from a NLWG has previously 
been reported as a route by which unstable nonlinear guid¬ 
ed waves (NLGW’s) decay. 10 However, in that case the 
input flux Fixes the input beam profile whereas here we 
hold the input proFile Fixed which is more in line with ex¬ 
perimental procedures. Potential applications of this ef¬ 
fect include optical limiting, coupling of adjacent 
waveguides in the absence of evanescent Field overlap, and 
a light-driven angular scanning element. 

We consider TE waves of frequency to in a slab 
waveguide, then assuming that the electric Field is homo¬ 
geneous in the y direction >z and x being the propagation 
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and transverse coordinates in units of c/to ) and writing 
the Field as 

E(r,/)= jy[ IV{x,z)e''0* ~‘'"' + c.c. ] (1) 

yields the usual slowly varying envelope equation for 

IPIX.Z):" 

2ip~ + — 1 % -[p : -n : ix, i W : )]W=0 . (2) 

dz a* - 

Here the refractive index in the various media is given by 

n 2 ix, W 2 ) = n y +a Y IF l x,z) 2 , y=c,f,s (3) 

the subscripts c.f.s referring to the cladding (,t< -dK 
Film {—d<x<d), and substrate ix>dK respectively. 
For the results presented here we specifically set 
n c =n j = 1.55, tty =1.57, a c = 10 : , and «y=a,--0, 
which corresponds to a self-focusing Kerr-type nonlinear 
cladding. 

The NLGW's are found as the solutions of Eq. (2) such 
that IVix.z)- R’dlx): 12 Equation (2) along with the con¬ 
tinuity conditions on the tangential components of the 
electric and magnetic Fields can then be used to generate 
the dispersion curves (guided wave flux S versus effective 
index P) for the system. In Fig. 1(a) we show the disper¬ 
sion curve for the TE,, NLGW's of a NLWG described by 
the parameters given above (note that since there is only 
one nonlinear medium the NLGW is uniquely specified 
by P). For this particular configuration the negative-slope 
region of the dispersion curve is unstable under propaga¬ 
tion and is indicated by a dashed line, the regions marked 
I and II are stable. 11 The inset beside each branch shows 
the general nature of the solution on that branch ie.g.. on 
branch II the solution is concentrated around the film¬ 
cladding interface 1 '), and the critical flux Y ( is that lluv 
above which the NLGW's on branch l cease to exist lot 
the example here S, " 0 195 i. 

Our numerical experiment consisted of launching an in 
pul beam of fixed profile but variable flux onto the 
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Some features of this soliton emission draw attention. 
Firstly and perhaps most importantly, the first threshold 
provides a highly desirable optical limiter characteristic: 
the input-output transfer is essentially 100% right up to 
this point at which the output drops drastically. Presum¬ 
ably the switch contrast shown in Fig. Kb) could be fur¬ 
ther increased by varying the guide parameters. Further¬ 
more, the output, which is always a branch-I NLGW and 
closely resembles the linear TE„ mode of the waveguide, is 
collimated. Thus, in contrast to previous proposals for 
optical limiters based on nonlinear refractive effects, 15-17 
the present system does not impose any undesirable 
phase-curvature on the output beam. Secondly, the angle 
of emission of the solitons is a function of the input flux. 
This effect was also observed by Tomlinson et at. in 
their study of nonlinear interfaces. As pointed out by 
these authors, this effect has potential use as a light- 
driven angular scanning element. Finally, we have found 
that a soliton emitted from one NLWG can be at least 
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partially trapped by a second adjacent guide. This offers 
the possiblity of directional-coupler action without the 
need for evanescent field overlap of the guides. Further¬ 
more, since soliton emission occurs only above a given 
threshold flux, coupling between the guides will experi¬ 
ence the same threshold behavior which could be used to 
construct logic elements. 

In conclusion, we have shown that external excitation 
of a NLWG can produce sequential threshold behavior 
through multisoliton emission, and a simple explanation 
for the threshold behavior has been advanced based on the 
NLGW dispersion curves. 

This research was supported by the Army Research Of¬ 
fice (ARO) (DAAG29-85-K-0026), the Air Force Office 
of Scientific Research (AFOSR), the Joint Services Optics 
Program of ARO and AFOSR, and the Science and En¬ 
gineering Research Council, U.K. (SERC GR/D84726). 


Opt. Lett. II. 315 (1986). 

"J. V. Moloney. J. Ariyasu, C. T. Seaton, and G. 1. Stegeman, 
Appl. Phys. Leu. 48. 826(1986). 
i: C. T. Seaton. J. D. Valera. R. L. Shoemaker. G. 1. Stegeman. 
J. Chilwell, and S. D. Smith. J. Quantum Electron. QE-21, 
774 (1985). 

K, H. Vach, C. T. Seaton, G. I. Stegeman, and I. C. Khoo, Opt. 
Lett. 9, 238(1984). 

l4 We have shown that in order to selectively excite the NLGW’s 
one must tailor the input beam to match that solution; failure 
to do so leads to decay of the input field: E. M. Wright, G. I. 
Stegeman, C. T. Seaton, and J. V. Moloney, Appl. Phys. Lett. 
49, 435 (1986). 

i? J. A. Hermann and M. E. Grigg, Opt. Commun. 49. 367 
(1984). 

ifl J. A. Hermann, J. Opt. Soc. Am. B I. 729 1 1984). 

,7 E. W. Van Stryland, H. Vanherzeele, M. A. Woodall, M. J. 
Soileau, A. L. Smirl. S. Guha. and T. F. Boggess. Opt. Eng. 
24, 613 (1985). 





























